pub fn new_birthday_probability(n: u32) -> f64 {
    let days = 365f64;
    let mut temp1 = 1f64;
    
    for i in 1..n {
        temp1 *= (days - i as f64)/days;
    }

    ((1f64 - temp1)* 10000f64).round() / 10000f64
}
